AD-A175  472 
UNCLASSIFIED 


GENERATION  OF  A  COMPOSITE  GRID  FOR  TURBINE  FLOUS 
CONSIDERATION  OF  A  N  CU>  NATIONAL  AERONAUTICS  AND 
SPACE  ADMINISTRATION  CLEVELAND  OH  LE  V  CHOO  ET  AL 
NOV  86  NASA-E-3201  NASA-TM-88890  F7G  20/4 


AD-A175  473 


NASA 

Technical  Memorandum  88890 


USAAVSCOM 
Technical  Report  86-C-38 


Generation  of  a  Composite  Grid  for 
Turbine  Flows  and  Consideration 
of  a  Numerical  Scheme 


Y.  Choo 

Lewis  Research  Center 
Cleveland,  Ohio 

and 


S.  Yoon 

Sverdrup  Technology,  Inc. 
Lewis  Research  Center 
Cleveland,  Ohio 

and 


C.  Reno 

Propulsion  Directorate 

U.  S.  Army  Aviation  Research  and  Technology  Activity— A  VSCOM 
q  Lewis  Research  Center 
CJ>  Cleveland,  Ohio 


LxJ 

. — I 

l  - 

,  C--3 

November  1986 

r  ' 


thia  iocoxnant  haw  ba«a 

ot  pobllc  rm*\  m 

(Jibuti on  la  riHrrJtad 


DTIC 

ELECTE 

DEC  3  0  1986 

I 

E 


GENERATION  OF  A  COMPOSITE  GRID  FOR  TURBINE  FLOWS  AND 


CONSIDERATION  OF  A  NUMERICAL  SCHEME 
Y.  Choo 

National  Aeronautics  and  Space  Administration 
Lewis  Research  Center 
Cleveland,  Ohio  44135 

S.  Yoon 

Sverdrup  Technology,  Inc. 

Lewis  Research  Center  /  01 
Cleveland,  Ohio  44135  (  c°pY 

rri 


Accession  For 

NTIS  GRA&I  ij/ 

DTIC  TAB 

Unannounced  q 

Justification _ 

By— _ 

Distribution/ 
Availability  Codes 
Avail  and/or 
list  Special 


C.  Reno* 

Propulsion  Directorate 

U.S.  Army  Aviation  Research  and  Technology  Activity  - 

Lewis  Research  Center 
Cleveland,  Ohio  44135 


SUMMARY 


AVSCOM 

Criminal  ocrtr -irvs  color 
F  ■  A*. I  ri’IC  roproduct- 

srill  be  in  black  and 
!te» 


A  composite  grid  was  generated  for  flows  in  turbines.  It  consisted  of  the 
C-grid  (or  0-grid)  in  the  immediate  vicinity  of  the  blade  and  the  H-grid  in  the 
middle  of  the  blade  passage  between  the  C-grids  and  in  the  upstream  region. 

This  new  composite  grid  provides  better  smoothness,  resolution,  and  orthogo 
nallty  than  any  single  grid  for  a  typical  turbine  blade  with  a  large  camber  and 
rounded  leading  and  trailing  edges.  The  C-H  (or  0-H)  composite  grid  has  an 
unusual  grid  point  that  is  connected  to  more  than  four  neighboring  nodes  in  two 
dimensions  (more  than  six  neighboring  nodes  in  three  dimensions).  A  finite- 
volume  lower-upper  (LU)  implicit  scheme  to  be  used  on  this  grid  poses  no 
problem  and  requires  no  special  treatment  because  each  interior  cell  of  this 
composite  grid  has  only  four  neighboring  cells  in  two  dimensions  (six  cells  in 
three  dimensions).  The  LU  Implicit  scheme  was  demonstrated  to  be  efficient  and 
robust  for  external  flows  in  a  broad  flow  regime  and  can  be  easily  applied  to 
Internal  flows  and  extended  from  two  to  three  dimensions. 


INTRODUCTION 

The  generation  of  a  good  grid  is  essential  to  obtaining  by  numerical  simu 
lation  an  accurate  solution  for  the  complex  fluid  flow  phenomena  in  a  turbine. 

A  smooth  and  nearly  orthogonal  H-grld  can  be  generated  if  the  turbine  blades 
are  thin  and  have  only  a  slight  camber  and  sharp  leading  and  trailing  edges. 

But  such  geometry  is  rare  for  a  turbine  because  of  the  high  work  factor  desired 
from  the  turbine  and  the  mechanical  factors  to  be  considered  in  designing  the 
blades,  which  usually  operate  at  high  temperature  and  speed  (ref.  1).  Turbine 
blades  are  often  designed  to  have  substantial  thickness  and  camber  and  rounded 
leading  and  trailing  edges.  Geometries  are  more  complicated  in  the  radial- 
inflow  turbine  than  in  the  axial  turbine.  A  typical  radial- inf  low  turbine 
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rotor,  for  Instance,  turns  the  flow  In  90°  deflection  In  the  meridional  plane 
while  simultaneously  turning  the  flow  In  the  blade-to-blade  plane.  Two  repre¬ 
sentative  turbine  geometries  are  used  to  examine  the  grid  generation  for  tur¬ 
bines:  a  radial-inflow  turbine  rotor  (e.g.,  ref.  2)  and  an  annular  turbine 
stator  (e.g.,  ref.  3). 

For  a  typical  turbine  blade  no  single  grid  offers  satisfactory  grid  prop¬ 
erties  In  the  entire  rotor  or  stator  passage.  With  the  H-grld  It  Is  difficult 
to  obtain  a  good  boundary-layer  resolution  around  the  leading  edge.  If  very 
fine  meshes  are  used  to  resolve  this,  the  number  of  grid  points  In  the  upstream 
region  becomes  excessive.  Thus  mesh  points  would  be  wasted.  With  the  standard 
C-grld  (or  0-grld),  for  a  typical  turbine  geometry,  the  grid  becomes  very 
skewed  because  the  blades  are  highly  cambered  and  the  blade-to-blade  grid  lines 
curve  sharply  In  the  midpassage  region  (refs.  4  and  5). 

A  two-dimensional  0/H  patched  grid  used  In  a  turbine  cascade  computation 
(ref.  6)  shows  the  slope  discontinuity  at  the  0-H  grid  Interface  and  has 
extremely  large-aspect-ratio  meshes  near  the  leading  edge.  The  numerical 
scheme  used  In  that  work  Is  a  cell-centered  scheme  based  on  the  Beam  and  Warm¬ 
ing  approximate  factorization  (ref.  7).  Although  the  flexibility  of  the  scheme 
Is  Indicated  In  two  dimensions,  Its  likely  limitations  In  three  dimensions  sug¬ 
gest  an  alternative  approach. 

In  this  study  two-  and  three-dimensional  composite  grids  were  generated 
In  an  attempt  to  Improve  the  grid  quality  In  terms  of  smoothness,  resolution, 
and  orthogonality.  A  numerical  scheme  that  will  run  on  this  composite  grid  Is 
discussed. 


STANDARD  GRIDS  FOR  TURBINE  GEOMETRIES 

A  three-dimensional  Invlscld  analysis  using  a  H-grld  was  performed  for  the 
radial-inflow  turbine  rotor  (ref.  8).  The  rotor  used  In  the  analysis  Is  shown 
In  figure  1.  The  rotor  had  12  full  blades  thick  enough  to  allow  Internal  cool¬ 
ing  passages  and  tralllng-edge  coolant  ejection.  The  rotor  had  unswept  radial 
blades  at  the  Inducer  Inlet  and  large  blockages  resulting  from  the  thick  trail¬ 
ing  edges.  In  figure  2  the  rotor  surface  Is  defined  by  using  some  of  the  grid 
lines  on  the  rotor  surface.  Three-dimensional  meshes  In  the  radial  turbine 
passage  are  shown  In  figure  3.  As  an  Illustration,  only  a  selected  number  of 
grid  lines  on  a  selected  number  of  grid  surfaces  are  shown  In  the  figure.  A 
detailed  H-grld  for  a  blade-to-blade  surface  Is  shown  In  figure  4(a).  The  mesh 
on  a  meridional  plane  Is  shown  In  figure  5(a).  Nonuniform  grid  spacing  was 
used  In  all  three  directions.  The  quasl-orthogonal  mesh  lines  were  clustered 
near  the  leading  and  trailing  edges,  as  Illustrated  In  figures  4(a)  and  5(a). 

In  addition,  this  feature  was  also  used  to  position  the  Initial  upstream  grid 
at  any  desired  location  with  a  selected  number  of  orthogonal  grid  lines  and  to 
position  the  final  downstream  grid  where  the  experimental  data  were  measured  so 
that  outputs  from  the  Invlscld  calculation  could  be  directly  compared  with 
experimental  results.  In  the  H-grld,  cusps  were  required  at  the  leading  and 
trailing  edges  to  get  convergence  of  flow  solution.  Figure  4(a)  shows  the  grid 
skewness  for  the  leading  edge;  figure  4(b)  shows  the  rapid  velocity  change  near 
the  leading  edge.  In  spite  of  the  difficult  blade  geometry  (l.e.,  thick  edges, 
large  trailing  blockage,  turning  of  the  passage  In  the  meridional  and  blade-to- 
blade  directions  simultaneously,  etc.),  the  three-dimensional  Invlscld  analysis 
predicted  the  flow  field  reasonably  well  In  a  qualitative  sense.  But  a  large 


entropy  change  occurred  around  the  leading-edge  region.  As  an  example  of 
results  obtained  from  the  three-dimensional  Invlscld  analysis,  figure  5(b) 
shows  secondary  velocity  vectors  on  cross-channel  surfaces. 

A  stator  vane  ring  Is  shown  In  figure  6(a).  The  full  annular  ring  con¬ 
sisted  of  36  vanes.  A  sector  of  four  vanes  used  to  make  laser  anemometer  meas¬ 
urements  Is  shown.  This  vane  ring  was  used  for  an  annular  cascade  experiment 
reported  In  reference  3.  The  stator  vane  geometry  Is  shown  In  figure  6(b). 

The  vanes  were  untwisted,  had  a  constant  profile  from  hub  to  tip,  and  had  a 
height  of  38.10  mm  and  an  axial  chord  of  38.23  mm.  The  vane  aspect  ratio  and 
the  solidity  at  the  mean  radius  (based  on  axial  chord)  were  1.0  and  0.93, 
respectively.  The  stator  hub-tip  radius  ratio  was  0.85  and  the  tip  diameter 
was  508  mm.  The  stacking  axis  of  the  vane  was  located  at  the  center  of  the 
tralllng-edge  circle.  To  alleviate  numerical  errors  associated  with  the  H-grld 
skewness  for  the  blunt  leading  edge,  a  standard  C-grld  was  generated  to  examine 
any  possible  advantage  over  the  H-grld  for  this  turbine  vane  geometry 
(fig.  7(a)).  The  grid  shows  large  skewing  on  the  suction  side  of  the  blade 
passage  (fig.  7(b)).  Results  of  numerical  simulations  using  this  grid  are 
likely  to  be  affected  by  the  grid  skewing,  as  Indicated  In  references  5  and  6. 


COMPOSITE  GRID 

Two-  and  three-dimensional  composite  grids  were  generated  for  the  turbine 
vanes  of  reference  3  (figs.  8(a)  and  (b),  respectively).  The  blades  were 
untwisted  and  stacked  at  the  trailing  edge  (fig.  8(c)).  Only  a  selected  number 
of  grid  lines  are  shown  for  Illustration.  This  composite  grid  consists  of  the 
C-grld  (or  0-grld)  In  the  Immediate  vicinity  of  the  blades  and  the  H-grld  In 
the  upstream  region  and  In  the  middle  of  the  blade  passage  between  the  C-grlds. 
The  C-grld  (or  0-grld)  portion  can  be  generated  by  using  either  the  elliptic 
method  (ref.  9)  or  the  algebraic  method  (ref.  10).  At  the  C-grld  and  H-grld 
Interface  (or  0-H  Interface)  the  slope  continuity  was  preserved  so  that  no 
special  numerical  approximations  are  needed  for  the  derivatives  at  the  Inter¬ 
face.  The  C-grld  (or  0-grld)  Is  orthogonal  to  the  blade  surface  and  provides  a 
good  boundary-layer  resolution  near  the  leading  edge.  This  composite  grid  has 
better  smoothness,  resolution,  and  orthogonality  than  any  type  of  single  grid 
for  a  typical  turbine  blade  with  large  camber  and  rounded  leading  and  trailing 
edges . 

For  the  two-dimensional  grid  In  figure  8(a)  the  C-grld  In  the  vicinity  of 
the  blade  was  generated  by  the  elliptic  grid  generation  code  (ref.  11).  Only  a 
portion  of  the  C-grld  near  the  blade  was  retained  for  the  composite  grid.  In 
the  C-grld  shown  In  figure  9,  for  Instance,  the  portion  between  n  =  0  and 
n  =  nc  <  nmax  was  retained,  where  n  »  0  and  n  *  n^ax  are  ^nner  and  outer 
boundaries,  respectively,  and  choice  of  nc  arbitrary  and  depends  on  the 
extent  of  the  shear  flow  region.  Then  the  H-grld  was  smoothly  patched  with  the 
C-grlds  on  the  pressure  and  suction  sides.  At  the  C-H  Interface  the  slope  con¬ 
tinuity  was  preserved.  The  cubic  spline  and  a  stretching  function  were  used  to 
generate  the  H-grld  In  the  midpassage  and  upstream  regions.  The  three- 
dimensional  grid  In  figure  8(b)  was  constructed  algebraically  from  the  two- 
dimensional  composite  grids  of  the  hub  and  shroud  surfaces. 


NUMERICAL  SCHEME 


Semidiscrete  Finite-Volume  Scheme 


The  two-dimensional  C-H  composite  grid  has  an  unusual  grid  point  that  Is 
connected  to  more  than  four  neighboring  nodes  (fig.  10(a)).  At  this  point  the 
usual  differencing  techniques  cannot  be  applied.  If  a  standard  finite- 
difference  scheme  were  used  to  solve  a  flow  problem  on  this  composite  grid,  It 
would  be  difficult  to  treat  this  special  point.  A  finite-volume  scheme 
(ref.  12)  to  be  used  with  this  composite  grid  presents  no  problem  and  requires 
no  special  treatment  because  each  Interior  cell  of  this  composite  grid  has  only 
four  neighboring  cells  In  two  dimensions  and  six  neighboring  cells  In  three 
dimensions  (fig.  10(b)).  The  finite-volume  scheme  Is  described  briefly  here. 


The  Euler  equations  In  Integral  form  can  be  written  as 


ft  fffv  *  Iff  •  is  ■  0 

n  an 


(l) 


for  a  fixed  region  n  with  boundary  an.  Here  w  represents  the 
conserved  quantity,  F  Is  the  corresponding  flux  term,  and  t  Is  time. 

A  convenient  way  to  ensure  a  steady-state  solution  Independent  of  the  time 
step  Is  to  separate  the  space  and  time  discretization  procedures.  In  the  semi¬ 
discrete  finite-volume  scheme  one  begins  by  applying  a  semidiscretization  In 
which  only  the  spatial  derivatives  are  approximated.  To  derive  a  semidiscrete 
model  that  can  be  used  to  treat  complex  geometric  domains,  the  computational 
domain  Is  divided  Into  quadrilateral  cells.  Assuming  that  the  dependent  vari¬ 
ables  are  known  at  the  center  of  each  cell,  a  system  of  ordinary  differential 
equations  Is  obtained  by  applying  equation  (1)  separately  to  each  cell.  These 
have  the  form 


ft  <Vu>  *  • 0  <2> 

where  S^j  Is  the  cell  area  and  Qij  Is  the  net  flux  out  of  the  cell. 

This  can  be  evaluated  as 

4 
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where  fk  and  g^  denote  values  of  the  flux  vectors  f  and  a  on  the  kth 
edge,  Axk  and  Ayk  are  the  Increments  of  x  and  y  along  the  edge  with 
appropriate  signs,  and  the  sum  Is  over  the  four  sides  of  the  cell.  The  flux 
vectors  are  evaluated  by  averaging  the  values  In  the  cells  on  either  side  of 
the  edge: 


fl  *  2  (fHlJ  +  f1J) 

for  example.  The  scheme  constructed  In  this  manner  reduces  to  a  central 
difference  scheme  on  a  Cartesian  grid  and  Is  second-order  accurate  In  space 
provided  that  the  mesh  Is  smooth  enough.  It  also  has  the  property  that  uniform 
flow  Is  an  exact  solution  of  the  difference  equations. 


LU  Implicit  Scheme 


The  composite  grid  Is  generated  to  be  used  for  both  Invlscld  and  viscous 
flow  calculations.  For  viscous  calculations  the  C-  or  0-meshes  must  be  very 
fine  to  resolve  the  boundary  layer.  And  It  Is  likely  that  the  time  step 
Imposed  by  an  explicit  stability  bound  will  be  much  less  than  that  Imposed  by 
the  accuracy  bound  of  an  Implicit  scheme.  Since  an  obvious  way  to  accelerate 
convergence  to  a  steady  state  Is  to  Increase  the  size  of  the  time  step,  an 
Implicit  scheme  Is  expected  to  have  a  faster  convergence.  Although  the  alter¬ 
nating  direction  Implicit  (AOI)  scheme  has  been  valuable  In  two-dimensional 
problems.  Its  Inherent  limitations  In  three  dimensions  suggest  an  alternative 
approach.  A  LU  Implicit  scheme  (ref.  13)  was  demonstrated  to  be  efficient  and 
robust  for  external  flows  In  a  broad  flow  regime  and  can  be  readily  extended 
from  two  to  three  dimensions.  This  scheme,  briefly  discussed  here,  will  be 
used  for  the  numerical  simulation  of  turbine  flows. 

The  conservation  law  form  of  the  Euler  equations  In  Cartesian  coordinates 
for  two-dimensional  flow  Is 


3W  3F  3G 
3t  *  3x  +  3y  *  8 

where  W  Is  the  vector  of  dependent  variables  and 
flux  vectors: 


(4) 

F  and  G  are  convective 


W  =  (p,pu,pv,pE)T 

2  T 

F  =  [pU.pU  +  p,pvu,u(pE  ♦  p)]  (5) 

2  T 

G  =  [pv ,puv ,pv  +  p,v(pE  +  p)] 


where  p,  u,  v,  E,  and  p  are  density,  velocity  components,  total  energy,  and 
pressure.  The  pressure  Is  obtained  from  the  equation  of  state 

P  -  Piy  -  d[e  -  J  (u2  ♦  v2)]  (6) 

where  y  Is  the  ratio  of  heat  capacity  at  constant  pressure  to  heat  capacity 
at  constant  volume.  These  equations  are  to  be  solved  for  a  steady-state 
dW/dt  *  0,  where  t  denotes  time. 

Let  the  Jacobian  matrices  be 


3F 
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and  let  the  correction  be 


<w  .  w"*1  -  h" 


where  n  denotes  the  time  level. 


The  linearized  Implicit  scheme  for  a  system  of  nonlinear  hyperbolic 
equations  such  as  the  Euler  equations  can  be  formulated  as 


(7) 


j^I  +  0  At  (OxA  ♦  OyB)J  AW  +  At  R  *  0 
where  I  Is  the  Identity  matrix  and  R  Is  the  residual 

E  =  DxF(Wn)  ♦  DyG(Wn) 

Here  Ox  and  Oy  are  central  difference  operators  that  approximate  a/ax 
and  a/ay. 


If  B  =  1/2  the  scheme  remains  second-order  accurate  In  time;  for  other 
values  of  B  the  time  accuracy  drops  to  first  order.  The  unfactored  Implicit 
scheme  (eq.  (7))  produces  a  large  block-banded  matrix,  which  Is  very  costly  to 
Invert  and  requires  huge  storaoe.  An  unconditionally  stable  Implicit  scheme 
that  has  error  terms  at  most  of  order  ( At ) 2  in  any  number  of  space  dimen¬ 
sions  can  be  derived  by  LU  factorization 


[ 
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*>][ 
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AW 


+  At  R  =  0 
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where  Dj  and  Oy  are  backward  difference  operators  and  D$  and  DJ 
are  forward  difference  operators.  The  reason  for  splitting  Is  to  ensure  the 
diagonal  dominance  of  lower  and  upper  factors  as  well  as  to  make  use  of  the 
built-in  Implicit  dissipation. 


Here  A+,  A-,  B+,  and  B-  are  constructed  so  that  the  eigenvalues  of 
matrices  are  nonnegative  and  those  of  matrices  are  nonpositive. 


A+  -  l  (A  ♦  rAI),  A'  -  ~  (A  -  rAI) 

B+  =  l  (B  ♦  rBI),  B"  =  \  (B  -  rBI) 


where 
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Here  Xa  and  Xb  represent  eigenvalues  of  Jacobian  matrices.  Equation  (8) 
can  be  Inverted  In  two  steps.  The  LU  Implicit  scheme  needs  the  Inversion  of 
sparse  triangular  matrices,  which  can  be  done  efficiently  without  using  large 
storage.  This  scheme  has  only  two  factors  In  three  dimensions. 


CONCLUDING  REMARKS 


Two-  and  three-dimensional  composite  grids  were  generated  In  an  attempt  to 
Improve  grid  quality  In  terms  of  smoothness,  resolution,  and  orthogonality  for 
an  annular  turbine  cascade.  This  composite  grid,  which  has  a  C-type  (or 
0-type)  grid  In  the  Immediate  vicinity  of  the  turbine  blade,  provides  a  good 
boundary-layer  resolution  around  the  leading-edge  region  for  viscous  calcula¬ 
tion,  has  orthogonality  at  the  blade  surface  and  slope  continuity  at  the  C-H 
(or  0-H)  Interface,  and  controls  mesh  distribution  In  the  upstream  region  with¬ 
out  using  excessive  grid  points.  This  composite  grid  eliminates  the  undesir¬ 
able  qualities  of  a  single  grid  when  generated  for  a  typical  turbine  geometry. 


A  finite-volume  lower-upper  (LU)  Implicit  scheme  can  be  used  In  solving 
for  the  turbine  flows  on  the  composite  grid.  This  grid  has  a  special  grid  node 


J 


that  Is  connected  to  more  than  four  neighboring  nodes  In  two  dimensions  and  to 
more  than  six  nodes  In  three  dimensions.  But  the  finite-volume  approach  poses 
no  problem  at  the  special  point  because  each  Interior  cell  has  only  four 
neighboring  cells  In  two  dimensions  and  only  six  cells  In  three  dimensions. 

The  LU  Implicit  scheme  was  proved  to  be  efficient  In  a  broad  flow  regime  and  Is 
expected  to  yield  accurate  solutions  on  the  Improved  composite  grid. 
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(a)  Grid  in  meridional  plane. 


(b)  Secondary  velocity  vectors  on  cross-channel  surfaces, 
figure  5.  -  Evolution  of  passage  vortex  (from  three-dimensional  Invlscid  computation  i. 


(a)  Ring  showing  cutout  tor  laser  window. 


V/Vcr  -  a  231 


i 


Radius, 


5.  Of  mm  - 


Axial 

centerline 


edge  circle 


Radius. 
0. 89  mm 


a-67° 


V/Vcr  •  a  778 


(b>  Geometry  at  mean  section.  (V/Vcr  •  critical  velocity  ratio.  I 
Figure  6.  -  Core  turbine  stator  vane 


Hi  Two  dimensions. 

Figure  &  -  Composite  grid  for  turbine  stator  vines  of  annular  cascade. 


lb)  Three  dimensions. 


ii)  Grid  nodes  to  be  used  in  finite-difference  scheme. 
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